//SR measurement of the fundamental length scales in the vortex state of YBa2Cu3O6.60 

J. E. Sonier, J. H. Brewer, R. F. Kiefl 
D. A. Bonn, S. R. Dunsiger, W. N. Hardy, Ruixing Liang, W. A. MacFarlane, R. I. Miller and T. M. Riseman * 
TRIUMF, Canadian Institute for Advanced Research and Department of Physics and Astronomy, University of British 

Columbia, Vancouver, British Columbia, Canada V6T IZl 

D. R. Noakes, C. E. Stronach and M. F. White Jr. 
Department of Physics, Virginia State University, Petersburg, Virginia 23806 
( ) 



The internal field distribution in the vortex state of 
YBa2Cu3O6.60 is shown to be a sensitive measure of both 
the magnetic penetration depth \ab and the vortex-core ra- 
dius po- The temperature dependence of po is found to be 
weaker than in the conventional superconductor NbSe2 and 
much weaker than theoretical predictions for an isolated vor- 
tex. The effective vortex-core radius decreases sharply with 
increasing H , whereas \ab{H) is found to be much stronger 
than in NbSe2. 



The magnetic penetration depth A and the coherence 
length ^ are the two fundamental length scales in a super- 
conductor. A is directly related to the superfluid density 
Ug whereas ^ is the length scale for spatial variations in 
the superconducting order parameter. In recent years, 
measurements of A have provided strong evidence for 
unconventional pairing in the high- Tc materials in that 
there are line nodes in the superconducting energy gap. 
In particular, in both the Meissner and vortex states, A 
has been observed to increase linearly with temperature 
T and magnetic field H ||,|| at low T in single crys- 
tals of YBa2Cu3 06.95. 

Much less is known about the behaviour of ^. In fact 
until now there has been no measurement of ^ deep in 
the superconducting state of a high-Tc superconductor. 
The magnitude of ^ near the phase boundary has been 
estimated from measurements of the upper critical field 
Hc2 using Ginzburg-Landau (GL) theory. For a type-II 
superconductor for H near i?c2, the coherence length is 
related to Hc2{T) in GL theory by. 

Reliable estimates of ^ in YBa2Cu307_5 made in this way 
are extremely difficult because the value of at T = is 
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SO large (i.e. ~100 T). Consequently, Hc2 measurements 
are limited to high temperatures above T/Tc ^ 0.85, so 
that the low T behaviour of ^(T) can only be determined 
by extrapolation. The situation is further complicated 
in the high-Tc compounds by strong thermal fluctuations 
over a sizeable region near Tc which result in broad transi- 
tions and poor estimates of Hc2{T). Furthermore, Eq. (|l|) 
may not be valid for an unconventional superconductor. 

It is desireable to have direct measurements of the co- 
herence length, which in the vortex state is related to the 
size of the vortex cores. In particular, for a conventional 
superconductor ^~poj where po is the vortex-core radius 
Q . In principle both STM and /LtSR can be used to char- 
acterize the size of vortex cores and thereby determine ^ 
directly deep in the superconducting state. Recent STM 
and muon spin rotation (/iSR) R| measurements on 
NbSe2 show that po decreases with increasing i7, as a re- 
sult of the increased interaction between vortices. Similar 
studies on high-T^ materials have not yet been performed 
with STM. Pinning of the vortices due to surface rough- 
ness and oxygen vacancies eliminates the long-range or- 
der in the vortex lattice and results in variations in the 
electronic structure of the cores. Furthermore, the ef- 
fect of the discontinuity in the quasiparticle excitation 
spectrum at the surface is still not understood. 

On the other hand, /iSR provides information on the 
vortex cores in the bulk of the sample. As explained in 
Ref. 0, in a /iSR experiment po is related to the high- 
field cutoff of the measured internal field distribution. 
The spectral weight at the cutoff grows as the density 
of the vortices increases. A well defined cutoff was ob- 
served in NbSe2 (Ref. [Q), where po is several times larger 
thanin YBa2Cu3 6.95 (Refs. In YBa2Cu3 6.95 no 

clear signal from the vortex cores was visible below 3 T. 
The temperature dependence of ^ (T) was investigated in 
an earlier /iSR study of YBa2Cu3 6.95 at higher mag- 
netic fields 1^. Unfortunately the signal-to- noise ratio 
was poor due to the infiuence of the large magnetic field 
on the positron orbits and timing resolution. In addition, 
demagnetization effects likely contributed significantly to 
the measured field distribution since the sample consisted 
of nineteen small crystals, and A and ^ were assumed in- 
dependent of H in the fitting procedure. 

In this Letter we present a /iSR study of the oxy- 
gen deficient high-Tc superconductor YBa2Cu3O6.60- 
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Compared to the optimally-oxygenated compound, 
YBa2Cu3O6.60 has a smaller carrier concentration in the 
Cu02 planes, a reduced Tc and i?c2, and a correspond- 
ingly larger i^. This allows us now to report the first 
detailed study of the fundamental length scales A and ^ 
in a high-Tc superconductor. 

We studied two different YBa2Cu3O6.60 samples with 
identical transition temperatures (59 K). The first sam- 
ple (SI) was obtained by deoxygenating the three-crystal 
mosaic of YBa2Cu3 6.95 used in Refs. The twin 

boundary spacing was on the order of 1 /im in the bulk. 
The second sample (S2) was grown from a separate batch 
and consisted of two large single crystals with a total a-h 
plane surface area of 30 mm^. S2 was mechanically de- 
twinned, such that the twin boundary density was about 
an order of magnitude smaller than in SI. Measurements 
were performed on field-cooled samples using the M15 
and M20 surface muon beamlines at TRIUMF and the 
same apparatus as that used in Refs. 1^,^,01. 

The experimental muon spin precession signal was fit 
assuming the local field due to the vortex lattice at any 
point in the a-h is given by , 

B{p)^B^{\-h^)Y, rrr^^, (2a) 



with, 



= 2^2^G2(1 + &4)[i_ 26(1 -6)2], 



(2b) 



where Bq is the average magnetic field, b — B/Bc2, 
£,ab is the GL coherence length and Ki{u) is a modified 
Bessel function. This analytical model of the field profile 
agrees extremely well with the exact numerical solutions 
of the GL equations |]l0| at low reduced fields b — whereas 
6 < 0.02 for our measurements. The term uKi{u) cuts 
off the summation thereby removing the logarithmic di- 
vergence of field at the vortex cores in the conventional 
London model. The cutoff is done in a way which pre- 
serves circular symmetry around the vortex cores. We 
note that the results herein are qualitatively similar to 
that obtained using a modified London model for B{p) 
with a Gaussian cutoff factor (as in Ref. jij ) and nearly 
identical to that using a Lorentzian cutoff — where only 
the latter is strictly valid at low fields. The internal 
field distribution n{B) was convoluted with a Gaussian of 
width a to account for vortex-lattice disorder and nuclear 
dipolar fields. 

The summation in Eq. (^ is taken over all reciprocal 
lattice vectors G of a triangular vortex lattice — the struc- 
ture which minimizes the free energy for a conventional 
superconductor. Infrared reflectance measurements of A 
in zero-field show a small anisotropy in the d-b plane 
(Aa/Ab = 1.3) which will stretch the triangular lattice, 
leading to elliptical cores in which ^b/£,a = ^a/ ^b- It was 
shown in Ref. Q through a scaling argument, that the 
corresponding magnetic field distribution is identical to 
the isotropic case. 
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FIG. 1. The Fourier transforms of tlie muon spin preces- 
sion signal in YBa2Cu3O6.60 after field cooling in a magnetic 
field H ~0.75 T down to T = 0.04 and 0.84 Tc (inset). The 
dashed curve is the Fourier transform of the simulated muon 
polarization function which best fits the data and the shaded 
region is the residual background signal. 

In a d2.2_j^2-wave superconductor the magnetic field 
distribution in the core region can be fourfold symmet- 
ric Jl2| , p^ and twofold symmetric with a-b anisotropy 
[Q. However, theoretical models for a (i^2_j,2-wave vor- 
tex core contain too many parameters to be useful in fit- 
ting experimental data. Modelling the core region with 
circular (or elliptical) symmetry should be sufficient to 
characterize the changes in core size with field and tem- 
perature. In general one expects the symmetry of the 
dx2^y2-waYe vortex cores to distort the lattice from tri- 
angular symmetry, but only at high fields where the in- 
tervortex spacing is small. So far no experiments have 
imaged the vortex lattice in YBa2Cu3O6.60- 

Figure 1 shows the real amplitude of the Fourier trans- 
form of the muon precession signal in YBa2Cu3O6.60 for 
T = 0.04 Tc and H = 0.75 T (solid curve) and of the 
simulated muon polarization function which best fits the 
data (dashed curve). Above the lineshape is nearly a 
perfect Gaussian with a width entirely due to the nuclear- 
dipolar fields. Well below Tc the internal field distribu- 
tion is very similar to that previously observed in NbSe2 , 
where the vortex lattice is known to be triangular. A 
small peak due to a residual background signal is also 
visible. At 0.84 Tc the field distribution is no longer 
asymmetric (Fig. 1 inset). We attribute this qualitative 
change in the lineshape to melting of the vortex lattice, 
ie: a transition from continuous 3D vortex lines to 2D 
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FIG. 2. The temperature dependence of \~^{T,H) in 
YBa2Cu3O6.60 (S2) for applied fields of 0.5 (solid circles), 
0.85 (open circles) and 1.25 T (solid triangles). Inset: The 
T-dependence of \ab at 0.5 T. The solid line shows the mi- 
crowave measurements of A\ab{T) = XabiT) — Aai,(1.25 K) in 
zero field [|l^ assuming our value Aai,(1.25 K) = 1762 A. 

vortex "pancakes" which are uncorrelated between 
planes. Such a transition has been observed in 
/iSR studies of highly anisotropic Bi2.i5Sri.85CaCu208+i5 
(BSCCO) [ID but does not occur in YBa2Cu3 06.95 in 
similar fields because of stronger interplane coupling. 
We estimate the variation of the crossover temperature 
r,„ with magnetic field for O.b < H < 1.5 T to be 
Tra=Tc - aH, where a = 17(1) K/T. 

Figure 2 shows the temperature dependence of X^^ (T) 
below Tm for three magnetic fields. As previously ob- 
served in YBa2Cu3 06.95 and Lai.85Sro.i5Cu04 [ p^ , 
there is a strong linear-T dependence at low temperature 
which is independent of H. The inset of Fig. 2 shows a 
comparison between the T-dependence of Xab{T) at 0.5 T 
and microwave cavity measurements |l7| of AXab{T) in 
zero magnetic field. The excellent agreement confirms 
that the fitting procedure which assumes a triangular 
vortex lattice, introduces at most only a small system- 
atic error in the absolute value of A. This is reasonable 
since it has been shown theoretically that including addi- 
tional terms in the free energy of the vortex state produce 
only minor changes in the internal field distribution [ p^ . 

As in Ref. § we define po to be the radius at which 
the supercurrent density Js(p) = V x B(p) reaches its 
maximum value. Js{p) profiles were generated from fits 
of the data to the field profile of Eq. (||) . Figure 3 shows 
the T-dependence oi pq{T) for the same fields as in Fig. 2. 
The solid lines are fits to the linear relation 
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FIG. 3. The T-dependence of po in YBa2Cu3O6.60 (S2) 
for applied fields of 0.5 (solid circles), 0.85 (open circles) 
and 1.25 T (solid triangles). Inset: The T-dependence 
of po(T)/po(0) in NbSe2 at 0.19 T (open circles) and in 
YBa2Cu3O6.60 at 0.5 T normalized to Tc = 59 K (solid circles) 
and T* = 250 K (open triangles). 

po(T)==po(0)[l +/3T/TJ, where ^-0.23 is essentially 
independent of field. The inset of Fig. 3 shows the T- 
dependence of po at 0.5 T normalized to po at T = 0. The 
linear term is much weaker than that found in NbSe2 at 
0.19 T ||l|l, where /3'-1.2. Better agreement (see Fig. 3) 
is obtained if T is normalized to T* = 250 K, the temper- 
ature below which a pseudogap opens in the spectrum of 

low-energy excitations in YBa2Cu3 06. 60 (Ref. 0). One 
possible interpretation is that the pairing amplitude is es- 
tablished prior to the onset of long range phase order at 
Tc |2|] . In both materials the size of the vortex core does 
not decrease as steeply with temperature as expected in 
theoretical predictions for a s-wave or (ij,2_j^2-wave 
p3[ superconductor. However, these theories pertain to 
a single isolated vortex and do not account for vortex lat- 
tice effects. Thermal fluctuations of the vortices about 
their average positions result in a premature truncation 
of the high-field tail in the p,SK lineshape — which results 
in an overestimate of po that increases with T. However, 
thermal fluctutations are expected to be most important 
at high magnetic fields |23) or near T^, and do not ac- 
count for the weaker T-dependence relative to NbSe2. 

In Fig. 4, Xab, Po and k extrapolated to T=0 are plot- 
ted as a function of H. The magnitude of Xab determined 
in SI is significantly lower than that in S2. The differ- 
ence is likely a result of vortex lattice distortions in SI 
due to twin boundary pinning. This introduces a system- 
atic uncertainty in the determination of Xab- The RMS 
deviation of the vortices from their positions in a perfect 
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FIG. 4. The field dependence of (a) Xab (b) po and (c) 
Kab = \ab/iab extrapolated to r = in YBa2Cu3O6.60- The 
data for SI (twinned) are shown as open circles whereas the 
solid circles designate S2 (detwinned). 

triangular lattice (estimated from the fitted values of 
(t) was found to be ~8% and 5% of the intervortex spac- 
ing for SI and S2, respectively. This disorder was in- 
dependent of H . The lines in Fig. 4(a) are fits to the 
linear relation, \ab{Q,H) = Aa6(0, 0)[1 -I- ^H/Hc2]- As- 
suming Hc2^70 T, Aafc(0,0) = 1586 A and 7 = 6.6 in SI, 
and Xab{0, 0) = 1699 A and 7 = 5.0 in S2. The increase 
in Xab with H is comparable to that observed recently 
in YBa2Cu3 06.95 ||] and is considerably stronger than 
that reported in NbSe2 where 7 = 1.6 This difference 
can be attributed to an enhancement in pair breaking 
caused by the applied field for an energy gap function 
with nodes — analogous to the nonlinear Meissner effect 

il- 

Figure 4(b) shows how po decreases with increasing 
magnetic field. Twin boundaries appear to have a neg- 
ligible effect on the core size since there is good agree- 
ment between SI and S2. A physical interpretation is 
that the increased interaction between vortices with field 
"squeezes" the vortex cores causing a reduction in po. 
The microscopic theory for a conventional superconduc- 



tor predicts such behaviour |2^. It is conceivable, the 
result may be attributed to quantum fluctuations at low 
temperatures which become important for large H [ p6| . 
However, quantum fluctuations are likely to have a small 
effect given the observed field dependence of po in NbSe2 
(Refs. H,01). The essential point is that quantum fluc- 
tuations are expected to be negligible in a conventional 
superconductor where ^ is large |26|| . It is important to 
realize that our measurements in YBa2Cu3O6.60 extend 
over a very narrow range of H/Hc2 relative to those for 
NbSe2. We were limited to this field range by the melting 
transition for higher H and the small amplitude of the 
high-field cutoff for lower H. Consequently, no reduction 
in the rate of change of po with H (expected for larger 
fields) was observed. A linear fit to the data yields an 
intercept po(0,0) = 89.5 A and slope -28.2 A/T. 

Figure 4(c) shows the iJ-dependence of k = Xab/^ab- 
A linear fit gives an intercept k{0, 0) = 17.0 and slope 
15.3 T^^. We also find that k is essentially indepen- 
dent of T over the entire temperature range. This agrees 
not only with magnetization measurements in the vortex 
state of BSCCO for T/Tc > 0.43 but also with our 
previous measurements in NbSe2. 

In conclusion, we have measured the T and H- 
dependences of Xab and po in YBa2Cu3O6.60 for H <^Hc2- 
We find that Xab and po both vary linearly with T at low 
temperatures. The T-dependence for po is considerably 
smaller than that found in NbSe2. Also, Xab increases 
while Po decreases with increasing magnetic field. We at- 
tribute the field dependence of po to the compression of 
the vortices due to vortex-vortex interactions. 
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